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The variational principle for the real homogeneous Monge- Ampere equation 
in two dimensions is shown to contain three arbitrary functions of four variables. 
There exist two different specializations of this variational principle where the 
Lagrangian is degenerate and furthermore contains an arbitrary function of 
two variables. The Hamiltonian formulation of these degenerate Lagrangian 
systems requires the use of Dirac's theory of constraints. As in the case of 
most completely integrable systems the constraints are second class and Dirac 
brackets directly yield the Hamiltonian operators. Thus the real homogeneous 
Monge- Ampere equation in two dimensions admits two classes of infinitely many 
Hamiltonian operators, namely a family of local, as well as another family non- 
local Hamiltonian operators and symplectic 2-forms which depend on arbitrary 
functions of two variables. The simplest nonlocal Hamiltonian operator cor- 
responds to the Kac-Moody algebra of vector fields and functions on the unit 
circle. Hamiltonian operators that belong to either class are compatible with 
each other but between classes there is only one compatible pair. In the case 
of real Monge- Ampere equations with constant right hand side this compatible 
pair is the only pair of Hamiltonian operators that survives. Then the com- 
plete integrability of all these real Monge- Ampere equations follows by Magri's 
theorem. Some of the remarkable properties we have obtained for the Hamilto- 
nian structure of the real homogeneous Monge- Ampere equation in two dimen- 
sions turn out to be generic to the real homogeneous Monge- Ampere equation 
and the geodesic flow for the complex homogeneous Monge-Ampere equation 
in arbitrary number of dimensions. Hence among all integrable nonlinear evo- 
lution equations in one space and one time dimensions, the real homogeneous 
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Monge- Ampere equation is distinguished as one that retains its character as an 
integrable system in multi-dimensions. 



1 Introduction 

Recently we had found that the real homogeneous Monge- Ampere equation 
(RHMA) in two dimensions, hereafter to be referred to as RHMA2, admits gen- 
eralized Hamiltonian structure. In field theory this type of structure arises for 
degenerate Lagrangian systems where we must use Dirac's theory of constraints 

to cast the problem into Hamiltonian form and finds its zenith in the theorem 
of Magri |3| for completely integrable bi-Hamiltonian systems. We refer to refs. 
[Q for an exposition of this subject. In this paper we shall show that RHMA 
and the geodesic flow for the complex homogeneous Monge-Ampere equation 
(CHMA) admit infinitely many Hamiltonian structures in arbitrary dimension. 
By Magri's theorem this result shows the complete integrability of RHMA and, 
in turn, RHMA provides the richest illustration of Magri's theorem. 

It will be useful to first list the various homogeneous Monge-Ampere equa- 
tions that we shall discuss in order to fix the notation. On a real manifold M 
of dimension n, RHMA„ is given by 

deiUap — Q a = 0, 1, n — 1 (1) 

where 

is the matrix of second derivatives. We shall take x*^ = t and write eq.(0) as a 
first order system of nonlinear evolution equations. Without loss of generality 
we may assume 

detUab^O, a~l,...,n — l (2) 

which is a statement of the nondegeracy of eq. (Q) . In 1 -I- 1 dimensions we have 
the simplest case of real Monge-Ampere equations 

UttUxx-UtJ^ = -K, 
K = ±1, 

which is hyperbolic, elliptic, or homogeneous (RHMA2) respectively. We shall 
show that there exist two types of Hamiltonian operators for RHMA2 , a family of 
local (^9|), as well as nonlocal ( pT| ) operators, both of which contain an arbitrary 
function of two variables. 
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We shall start with a discussion of the variational principles underlying 
RHMA2. The Lagrangian which results in the equations of motion for RHMA2 
contains three arbitrary functions of four variables. We shall concentrate on 
two qualitatively different classes of Lagrangians which are specializations of 
the master Lagrangian (|lO|) for RHMA2, c.f. eqs.(^l]) and ( p3[ ) below, both of 
which contain arbitrary functions and which are furthermore degenerate. That 
is, the passage to a Hamiltonian formulation of these degenerate Lagrangian 
systems gives rise to second class constraints as in the case of most completely 
integrable systems The resulting two families of Hamiltonian operators, c.f. 



eqs.(|29D and (51) contain all the information in the Dirac brackets of RHMA2. 
These Hamiltonian operators are infinite in number as they contain arbitrary 
functions which make their first appearance in the Lagrangian formulation of 
RHMA2 which, in turn, is related to its character as a universal field equation |^ . 
Any pair of local, or nonlocal Hamiltonian operators are compatible, however a 
local Hamiltonian operator is not in general compatible with a nonlocal Hamilto- 
nian operator with the exception of Jq of eq.(|30|) and Ji given by eq. (|5^) . These 
two Hamiltonian operators determine part of the multi-Hamiltonian structure 
of the real Monge- Ampere equation (||) when the right hand side is a non-zero 
constant. 

We shall show that the Kac-Moody algebra that corresponds to the simplest 
nonlocal Hamiltonian operator (^2|) for RHMA2 consists of the algebra of vector 
fields and functions on S^. The family of local Hamiltonian operators can be 
brought to the standard form of a first order operator by a Miura transformation. 
The recursion operator obtained by a composition of the simplest local and 
nonlocal Hamiltonian operators is Sheftel'-type [Q, that is, it can be written as 
the square of a first order operator which is a good recursion operator itself. 
Finally, we shall present the symplectic structure of RHMA2 which is dual 
to these Hamiltonian structures. We shall also show that the symplectic 2- 
forms obtained in this way are the time components of the Witten-Zuckerman 
[D 2-form that is the main geometrical object in the covariant formulation of 
symplectic structure. 

In the theory of developable surfaces there is an equation which may be 
called Ur-RHMA2 

UtUx ^ k (4) 

where k is constant, as it is well-known |^ that solutions of cq.(^ must satisfy 
RHMA2 . We shall discuss the Hamiltonian structure of Ur-RHMA2 in section 
^ and show that the Hamiltonian operators appropriate to this equation consist 
of a scalar version of the Hamiltonian operators for RHMA2 which lends further 
credence to its name. 

We shall show that RHMA is an example of a nonlinear evolution equation 
that is integrable in arbitrary number of dimensions, see ^ for other universal 
field equations with this property. In section ^ we shall find that some of 
the remarkable results we shall report on the Hamiltonian structure of RHMA2 
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hold for RHMA„ as well. Hence RHMA2 will serve as the prototype for the 
discussion of the symplectic structure of all real homogeneous Monge- Ampere 
equations. The Hamiltonian structure of the complex Monge- Ampere equation 
is quite different from that of RHMA and will not be considered in this paper. 
On the other hand the Hamiltonian structure of the geodesic flow defined by 
CHMA which will be discussed in section is very similar to that of 
RHMA. 



2 First order evolutionary form of RMA2 

For the Hamiltonian treatment of real Monge- Ampere equations we need to 
rewrite them as a system of nonlinear evolution equations which is first order 
in time. For RHMA2 this is most conveniently accomplished by introducing the 
definitions 

Ux =p, ut^q (5) 
and accordingly the equations of motion are either given by 

Ut =q, qt = — [qj - K) , (6) 

Uxx 

or ^ 

Pt = qx, qt = — {qx^ - K) (7) 

Px 

but we note that this split of eq.(|3|) into either one of these pair of evolution 
equations is not unique. Thus we have necessarily introduced a degree of free- 
dom over and above that present in eq.(0) itself. We shall return to this point 



in section 5.S on the Miura transformation. 

Henceforth we shall use the notation {m*; i = 1,2} for the first order variables 
where will stand for either m, or p and v? = q. Then the vector field defining 
the flow for eq.(||) is given by 



-^u — q^ -\ -^p — qx^ ^ (o) 

du oq op Px oq 



respectively. We shall cast the equations of motion following from eqs.(g) into 
the form of Hamilton's equations 

u\ = X(u') ^ {u\ H} ^ J'^ SkH (9) 

where the Hamiltonian operator J deflning the Poisson bracket is a skew- 
symmetric matrix of differential operators satisfying the Jacobi identities and 
di denotes the variational derivative with respect to u'. 
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3 Variational formulation of RHMA2 



We shall find that the Hamiltonian operators for RHMA2 are doubly infinite 
in number. The reason for the existence of such a multitude of Hamiltonian 
structures can be traced back to the fact that the variational principle (^o|) 
for the real homogeneous Monge- Ampere equation contains arbitrary functions. 
Our approach to the construction of the Hamiltonian operators as Dirac brackets 
appropriate to degenerate Lagrangian systems will make it manifest that the 
existence of arbitrary functions in Lagrangians for RHMA2 is responsible for 
their appearance in the Hamiltonian operators. 

It can be directly verified that the equations of motion following from the 
variational principle 

SI = 0, I = ( C dtdx 



with Lagrangian density 

C ^ li Utt + h utx + h , 

/ \ (10) 

/ — I 1 11 11 ^tt Utx \ n — ^ 'I 
Iq, — \ utx^ Uxx J u; — -L,^,0 

yield RHMA2 in the form of eq.(||). There are three arbitrary functions of 
four variables in this Lagrangian. We have not explicitly indicated the possible 
dependence of la on the third combination of the ratio of second derivatives 
above as it is derivable from the others. It appears that this is the richest 
example of a variational principle for any nonlinear partial differential equation. 
Fairlie, Govaerts and Morozov have considered the hierarchy of field equations 
where the equations of motion at any level are proportional to the Lagrangian 
at the next. They pointed out that this is a finite hierarchy ending at "universal 
field equations." RHMA2 is precisely such a universal field equation and this is 
the reason why we should expect arbitrary functions in the Lagrangian (pO|). 

In the variational formulation of eqs.(p|we shall consider various specializa- 
tions of the arbitrary functions la in eq.(|lC|) which will still contain arbitrary 
functions of a more specialized variety. It is, however, important to remember 
that the enormous number of possibilities offered by the master Lagrangian ( p^ 
is the source of all the results we shall present below. 

4 Local Dirac bracket for RHMA2 

A useful specialization of the Lagrangian density ( p^ that results in the first 
order equations of motion (||) for RHMA2 is given by 

Cx^ Xqt- Xu^ [ut - q) qx (H) 

where 

X^Xiux,q), A„^,^^0 (12) 
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is a twice differentiable arbitrary function of two variables. The canonical mo- 
menta appropriate to this Lagrangian are given by 



71'2 



dut 



-X. 



dqt 



A _ 



= A 



subject to the canonical Poisson bracket relations 

{^,{x),u\y)} = 5l5{x-y) 



(13) 



(14) 



with all others vanishing. But the momenta (13) cannot be inverted for the 
velocities, or alternatively the Hessian vanishes 



det 



= 



(15) 



and we have a degenerate Lagrangian system. Thus the passage to the Hamil- 
tonian formulation of the Lagrangian (11) requires the use of Dirac's theory of 
constraints. 

Following Dirac we introduce the two primary constraints that result from 
eqs.(|l|) 

01 = TTm + Xu^ q-x , 

(16) 

(f>2 = TTq - X 

and calculate their Poisson brackets using the canonical Poisson bracket rela- 
tions (|l|). The resuh 

Wi{x),<l>i{y)} = Qy Xuyuy Sy{x - y) - Xu^u^ S^{y - x) , 

{(f>i{x),(j>2{y)} = Xuy Sy{x - y) - Xu^ 6x{y - x) - qx Xqu^ S{y - x) , (17) 

{(t>2{x),My)} = 

shows that the constraints (^6|) are second class as they do not vanish modulo 
the constraints. This is a typical situation for integrable systems as in the 
example of KdV, or shallow water equations ||]. I am grateful to C. A. P. 
Galvao [0 for pointing out to me that one should not simplify the Poisson 
brackets of the constraints (17) using the rules for manipulating distributions 
as such "simplifications" often lead to incorrect Dirac brackets. 
The total Hamiltonian of Dirac is given by 



Ht = J {^i u\- C + c"<j)i) dx 



(18) 



where c* are Lagrange multipliers and summation over i = 1, 2 is implied. The 
condition that the constraints are maintained in time 



{(b,{x),HT}^0 



(19) 
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gives rise to no further constraints which would have been secondary constraints. 
Instead, using eqs.(p^, we find that the Lagrange multiphers are determined 
from eqs.(pj|) 

ci = g, 0^ = ^ 

and they do not depend on the choice of the arbitrary function A. This is 
expected since the constraints and therefore the total Hamiltonian is linear in 
the momenta, the correct equations of motion will result only if the Lagrange 
multipliers are simply the components of the vector field (||) for the flow. Finally 
from eq.(p^ we flnd 

HxT^ f (q vTu + ^7:,- ^) dx (20) 

J \ ^xx ^xx J 

for Dirac's total Hamiltonian. The check that with this total Hamiltonian all 
the equations of motion are satisfled is straight-forward and we can summarize 
all of them in Hamilton's equations 

At = {AH)^T} (21) 

where A is any functional of the canonical variables. There is, however, one 
further simplification that we can carry out because in Dirac's theory second 
class constraints hold as strong equations. This fact gives us the choice of 
eliminating the momenta from eqs.(|l^). Then we can write 

UxT^qqxK^ (22) 

for the total Hamiltonian density of Dirac. 

Given any two differentiable functionals of the canonical variables A and 
the Dirac bracket is defined by 

{A{x),B{v)}d^ {A{x),B{y)} 

f (23) 

- / {A{x),M^)}J'\z, w){Mw),Biy)} dz dw 

where J^^ is the inverse of the matrix of Poisson brackets of the constraints. 
The definition of the inverse is simply 

{0,(x), 0fe(z)} J^^ (z, y) dz = Si 5{x - y) (24) 

which results in a set of differential equations for the entries of . Starting with 
the Poisson bracket relations ([l^) we find that eqs.(^) can be solved readily to 
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yield 

Jl\x,y)^ 0, 

'^UxUx '-''XX 

4H^,y)- X \ 5i^-y)^ ^^^^ 

^UxUx ^XX 

Jf{^,v)= ~ , r "^^^ 6{x-y) 



for the inverse of ( |l7|) . It will be convenient to rename the arbitrary function A 
as /i according to 

K^u^ = -tJ'q (26) 
for ease of calculations that will follow. 



5 Local Hamiltonian operators for RHMA2 

The transition from the Dirac bracket to the Hamiltonian operator is given by 

{u^{x),u\y)}D = J'\x, y) = .r\x)5{x - y) (27) 

and from eqs.(^^ and ( p5| ) it follows that the Hamiltonian operator correspond- 
ing to the degenerate Lagrangian (ni]) is simply 



J, 



g \ 

(J'q "^xx 

q qqx p -l n -Ms: 

II — j^x-r-Ux 2" / 



(28) 



which contains the arbitrary function /i of two variables. Two particular choices 
of this arbitrary function, namely /i = ^q^ and /i = Uxq result in the bi- 
Hamiltonian structure of eq.(||) reported in 

Under the change of variable p = Ux the Hamiltonian operator (|2^) is trans- 
formed 

( ' ""'^ ] 

which is the form of the local Hamiltonian operator for RHMA2 that we shall use 
henceforth. The most important element of the family of Hamiltonian operators 
in eq.(^9|) is given by 



■h - 



/ Dxj-^ 
\ -^Dx ^^Dx + Dx-^ 

\ Px Px 



(30) 
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which corresponds to the choice fi = . 

From the construction of the Dirac bracket in section ^ it is clear that the 
reason for the presence of the arbitrary function /i in the Haniihonian operator 
( p9[ ) can be traced back to the degenerate Lagrangian (|l^) where it makes its 
first appearance as A. In order to proceed we need to prove that (|2^) is indeed a 
Hamihonian operator which requires a check of the Jacobi identities. However, 
this also follows from the fact that (|2^) is derived from the Dirac bracket ( p5| ) 
for which we have a general proof of the Jacobi identities ||l^ . 



5.1 Jacobi identities 



The Hamiltonian operator is a bi- vector which defines the Poisson bracket. The 
skew-symmetry of the operator ( |2^ ) is manifest. In order to verify that a given 
bi-vector is Hamiltonian we must verify that it satisfies the tri-vector Jacobi 
identities. Thus following Olver Q we introduce an arbitrary basis of tangent 
vectors which are then conveniently manipulated according to the rules of 
exterior calculus. The Jacobi identities are given by the compact expression 



L SI — {mod. div. 



where 



L = je, 



I^-e'^ AJQ 
2 



(31) 



(32) 



and S denotes the variational derivative. The vanishing of the tri-vector ( pl| ) 
modulo a divergence is equivalent to the satisfaction of the Jacobi identities. 

For the Hamiltonian operator ( p9| ) we have a two component system and 
introducing the basis 

' V 



e 



we have from eqs.(p3) 



L,. 



\ IJ'qPx 



I o qqx Q , f qqx 



n^x 



(33) 



NPx Kf^qPxJ X / 
Qx 



'AVx 



l-^qPx \ Px 

To form the required expression in cq. (|3l|) we first calculate the variational 



9 



derivatives 



dp 



f^qPx h'-q Px H-qPx 



ri Qx + mxx o f^qq^Qx o^pgOT^ g qQxPxx 
Z Z 2 — ^ O 2 — 2 O 



' /\Vx 



mpq_ 

f^qPx 



) QPxx 
^NPx 



_q_ 



2 f\ Oxx 



NPx 

and take their exterior product with L^. The result 



Lfj. SI^ = - 



2 OAOx/\r]x 



(34) 



is a total derivative so that the Jacobi identities are satisfied. 



5.2 Hamilton's equations 

It can be directly verified that 

T-i^L = MPx (35) 

is conserved for the flow (||) and this is the Hamiltonian density appropriate to 
the operator (|29|). That is, the equations of motion (||) are cast into the form of 
Hamilton's equations with and H^. But since /i is an arbitrary function 
of p and q, given any other differentiable function ly = vij), q) we have 

5H^ = J, 5H, (36) 

which is a statement of the Lenard-Magri recursion relation. Thus we have 
another expression of the fact that there exist infinitely many Hamiltonian op- 



erators for RHMA2. The Casimir density for the Hamiltonian operator (29) is 
given by 

C^=crpx, (37) 
where ^ 

Oq = -fiq (38) 

q 

since it can be readily verified that J^5C^ = 0. The proof of the Jacobi identities 



for ( |29| ) in section 5.1 is for arbitrary /i, hence all local Hamiltonian operators 
are compatible with each other. 
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5.3 Miura transformation 



The remarkable feature of the Hamihonian operator ( pOj ) is the existence of /x, 
an arbitrary function of two variables. We have already remarked that for the 
Hamiltonian formulation of RHMA2 we need to start with a pair of nonlinear 
evolution equations and therefore had to use the equations of motion defined by 
the vector field (^) rather than eq.(||) itself. This introduces an extra degree of 
freedom as the definitions of p, q cannot be unique. It may appear that this extra 
degree of freedom could be responsible for the existence of the arbitrary function 
H in the Hamiltonian operator (|29|). This is not the case. We have shown 
that the reason for the appearance of this arbitrary function in the Hamiltonian 
operator can be traced back to the existence of arbitrary functions in the master 
Lagrangian. Nevertheless, there is certainly an extra degree of freedom available 
in the definition of the first order dependent variables for RHMA2 and it is 
possible to exploit it for various purposes. In fact by a redefinition of q the 
Hamiltonian operator ( p9| ) can be brought to the standard form of a first order 
operator with constant coefficients. 

In place of eq.(Q) we may use the following definitions for the auxiliary 
variables p, q 

Ux^P, ut = Q [p, — ] (39) 

V PxJ 

where Q is a differentiable function of its arguments. The equations of motion 
are now given by 

Pt = Qx, qt= ( — Qx) (40) 

\Px ) X 

which are already in standard Hamiltonian form. That is, we have 



Pt A / Dx\(i^ 
qt J \ Dx J\ 6, 



H (41) 



with the familiar first order Hamiltonian operator and the Hamiltonian density 
is given by 

n^h{p, j-^ Px (42) 

where h is related to Q through Hq = Q. The transformation of variables ob- 
tained by a comparison of eqs.(^) and (||) is a Miura transformation as it brings 
the Hamiltonian operator to the canonical form ( pl| ) of Gardner-Zakharov- 
Fadeev lIT^ with constant coefficients. 



6 Nonlocal Dirac bracket for RHMA2 

There is another family of Hamiltonian operators for eqs.(^ which falls outside 
the class of Hamiltonian operators in eq.(p9|). Its origin can be traced back to 
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the fact that 

-Cfc = HPx + Kr {pt - qx) (43) 
with K, an arbitrary function of two variables 

K = K{p,r) , r = — , Krr ^0 (44) 
Px 

is a new Lagrangian for the system (|^). This is another speciahzation of the 
master Lagrangian (|lO|) which is also degenerate and the passage to its Hamil- 
tonian formulation again requires the use of Dirac's theory. 

Starting with the Lagrangian (|4^) the momenta are given by 

ir 

Ha = = ^ = 



which cannot be inverted for the velocities. So we introduce the primary con- 
straints 

<i>i = n„ — Kr 

and using the canonical Poisson bracket relations (^4|) we find 

{$i(a;), $i(y)} Krrix) ^ Sx{x ~ y) - Krr{y) % Sy{y ~ x) 

Px Py 

{$i(a;),<i>2(y)} = ^rrix) ^ Sy{x - y) (47) 

{<^2ix),^2{y)} = 

which show that the constraints are once again second class. The total Hamil- 
tonian consists of a sum of the free Hamiltonian with a linear combination of 
the constraints as in eq.(^8|) and the requirement that the constraints are 
maintained in time, c.f. eqs.(p^), does not lead to secondary constraints but 
instead determines the Lagrange multipliers 

c'=qx, c' = ^ 

Px 

which also do not depend on the arbitrary function n. Using this information 
we find that ^ 

H,T = J (^qxttp + ^Uq- Kp^^ dx (48) 



is the total Hamiltonian and the equations of motion are given by eq.(21). Once 
again we can simplify the total Hamiltonian using the fact that second-class 
constraints (^) hold as strong equations and we find 

H^T ^ Krqx - HPx (49) 
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for the Hamiltonian density of Dirac. 

For the Dirac bracket we need the inverse of the Poisson brackets of the 
constraints (|4^). The solution of eqs.(]2^) is given by 

Jl\x,y)^ 0, 

7l2/^ „\ Px 



Jl\x,y)= -~5{x-y)j ^^di, (^0) 

where is the Heaviside unit step function. The Dirac bracket for the La- 
grangian (Esh now follows directly from eq.d23). 



7 Nonlocal Hamiltonian operators for RHMA^ 



The correspondence (|2^) enables us to express the result ( |50| ) for the Dirac 
bracket of the degenerate Lagrangian ( ^3|) in the form of the Hamiltonian oper- 
ator 

where is the inverse of Da;. We refer to jl^ for the definition and properties 
of D~^, in particular, 

and the integrals are taken in the principal value sense. This Hamiltonian 
operator is therefore nonlocal. Once again, the appearance of the arbitrary 
function n in the Hamiltonian operator ( |5l| ) is a consequence of the existence 
of arbitrary functions in the master Lagrangian (p^). 

The most important member of the family of nonlocal Hamiltonian operators 
is given by 

Ji = (52) 
which has a linear dependence on px and qx . It results from the Lagrangian 

A = (53) 
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with the simplest choice of arbitrary function k = ^r^. It will also be useful to 
rewrite the operator ( |52| ) for the system of variables consisting of u and q 

J. J " ] (54) 

which follows from the change of variable p — Ux- 



In section 7.3 we shall find that the nonlocal Hamiltonian operator ( p2[ ) 
has a natural interpretation, namely the Kac-Moody algebra corresponding to 
this operator is the algebra of vector fields and functions on the unit circle. 
This operator is also distinguished as the only element of the family of nonlo- 
cal Hamiltonian operators that survives in the case of the real Monge- Ampere 
equation with constant right hand side. 

7.1 Jacobi identities 

The proof of the Jacobi identities for of eq.(|5l|) proceeds along the gen- 



eral lines indicated in section 5.1 but the properties of Dx ^ must be carefully 



considered. For ease of writing we let p = (K^r) ^ and from eq.(pl|) we have 

/ PPxDx'^0 \ 

L.= \ , (55) 

V Dx-' (ppxv) + pq.Dx-'e + Dx-' (pqxO) J 

Ik = PPxV A Dx-'e + pqx9 A Dx-'e . 
The variational derivatives of are given by 



61. 
Sp 



(56) 



{rpr - p) {rix A D-'e + rj A 9) + rprxV A 
+ (rVr)^ 9 A Dx-'9 - r^pr9x A Dx-'9 , 

^ = ~PrxVADx-'9- pr{vxADx-'9 + rjA9) 

-{p + rpr) 9x A D -'9 - (p + rpr)x 9 A D-'9 

which results in a total derivative for the tri-linear form ( ^T| ) 

LJI. = - {Dx-\pxV + qx9) A [pii +{p + rpr) 9] A D-'9]^ 

completing the proof of the Jacobi identities. 

7.2 Hamilton's equations 

The Hamiltonian function i/^ appropriate to the nonlocal Hamiltonian operator 
yj) is given by a local expression which is simply the integral of the total 
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Hamiltonian density of Dirac in cq.(^) and once again we have the Lenard- 
Magri recursion relation 

J. SH, = J, SH, (57) 

where t is another arbitrary function of p, r. All nonlocal Hamiltonian operators 
are compatible with each other. 



7.3 Kac-Moody algebra 

Hamiltonian operators associated with integrable nonlinear evolution equations 
give rise to Kac-Moody (KM) algebras |]l5|. Two compatible Hamiltonian op- 
erators actually yield an infinite hierarchy of KM algebras. There is an explicit 
algorithm for the construction of KM algebras from the Hamiltonian operator 
which is essentially based on Fourier analysis. Since the nonlocal operator ( ^2[ ) 
depends linearly on and Fourier analysis makes sense and the operator 
( p2| ) is suitable for consideration as the backbone of a possible KM algebra. The 
nonlinearities in the local operator present a formidable obstacle to any similar 
discussion of (p9|). Using 

P{^) I -T'ne^^dn q{x) = / lQ„e*™-dm 

2n J n 2n J m 

we find that the KM algebra appropriate to the operator (|5^) is given by 

[^m ; ^n] — , 

[Pm^Qn] = mVm+n. (58) 
[Qm, Qn] = {m~n) Qm+n , 

which can be recognized as the algebra of functions and vector fields on 5*^. 
Thus we can use the representation 

V — 2^" O — 2^™+^ — 

dz 

for the algebra underlying the simplest nonlocal Hamiltonian operator (^2|) . 



7.4 Compatibility 

Two Hamiltonian operators are compatible if their linear combination with con- 
stant coefficients is also a Hamiltonian operator. In the case where one of the 
Hamiltonian operators belongs to the class of eq.(^9|) and the other one to 
eg. (|5l|) , the check of compatibility requires that 

C^,^ L^5h + 61^ = (59) 

modulo a divergence. It turns out that only the simplest local and nonlocal 



operators Jq and Ji are compatible. From eqs.(|33|), (p4), (55) and Ipu) it follows 
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that 




which estabUshes the compatibihty of Jq and Ji. In aU other cases local and 
nonlocal Hamiltonian operators are incompatible. 



8 Recursion operator 



Since there are infinitely many compatible Hamiltonian operators of both local 
( p9| ) as well as nonlocal ( |5l| ) variety, there exists infinitely many opportunities 
for constructing recursion operators of either type or both. In the first category 
we have 





7? — 7 7 — 



(60) 



which generalizes our earlier result while in the latter category we find 





7^K 



jj —1 2^. I ^rr 

^ L Px \ ^rr 

These recursion operators satisfy the Lax equation 



r) — 1 l^rr r) 



Tit = [7^, ^] 



where 



A = 



9: 



Px 



Bx 

o2x_ JJ 
^ n ^x 
Px 



(61) 



(62) 



(63) 



is obtained from the Frechet derivative of the flow (||). 

However, there is also the possibility of constructing recursion operators 
by the composition of Hamiltonian operators of both local and nonlocal types 
provided they are compatible. We have found that Jq and J\ form a compatible 
pair. Thus we are led to the recursion operator 



I 



\ 



_ J_ 

'^Px^'^Px 



Px Px 

and we find that a remarkable factorization takes place 







Px^'-'Px^'' 



(64) 



(65) 
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where 







^=1 1 (66) 



is a first order operator. Tliis situation is familiar from third order Hamiltonian 
operators for two component equations of hydrodynamic type . The recursion 
operator obtained from a composition of third and first order Hamiltonian op- 
erators factorizes to yield the first order recursion operator of Shcftcl'. Similar 
considerations suggest that £ itself is a good recursion operator. It is readily 
verified that £ as well as its inverse 

£-^=( ' ] (67) 

satisfy eq.(|62|) and are therefore good recursion operators. Finally, we should 
remark that the recursion operators ( |60|) and (^Tj) can also be factorized, how- 
ever, unlike the case of eq.(|65|), the factors consist of different operators and 
each one is not a recursion operator. 

8.1 Infinite sequences of Hamiltonians 

There are several infinite families of conserved quantities that one can obtain 
by repeated application of the recursion operators we have presented in section 
Ij. Corresponding to the recursion operator £ we have J- which is defined by 

S^J„-^Ji=TT (68) 

which yields infinite sequences of gradients of conserved quantities. We find 
that J-" is given by 



= ^ ^ (69) 




and a simple example of the infinite sequence of conserved quantities that it 
generates is the following 



... ^ Ip-g,. (70) 



n 



Depending on the starting point one can construct infinitely many such se- 
quences. For example, S generates a family of local as well as nonlocal conserved 
quantities according to the scheme 



S 



{^^^ [^^^ (pf)]} ^ ^ipt 

(71) 

^^q^Px^ ^qPxDjT^ [pxD^^ (qp^)] 
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and there are infinitely many different combinations of such sequences. The 
above expressions for the conserved Hamiltonian densities can be simphfied 
by discarding divergences, however, this form is useful because it suggests the 
generic term for the conserved quantity obtained by repeated applications of 
<S. The determination of the cardinality of the conserved Hamiltonians which 
are in involution with respect to Poisson brackets defined by both types of 
Hamiltonian operators is complicated because of the appearance of arbitrary 
functions in these Hamiltonian operators. 

8.2 Higher flows 

Starting with the flow (H) we should obtain higher flows for RHMA2 through 
the application of any one of the recursion operators (|60|), (|6l|), (|66|), or ( |67|) 
on the vector fleld (||). However, in all of these cases an elementary calculation 
shows that the resulting higher flow is simply 

PtQx - qtPx ^ , 
the real homogeneous Monge- Ampere equation back again. 

9 Symplectic form of RHMA2 

The principal geometrical object in the theory of symplectic structure is the 
symplectic 2-form uj which is closed 

duj = 0, (72) 

and by Poincare's lemma lo can be written as 

u! = da (73) 

in a local neighborhood. On the other hand the Hamiltonian operator maps 
differentials of functions into vector flelds which works in the opposite direc- 
tion. Thus the statement of the symplectic structure of the equations of motion 
consists of 

ixuj = dH (74) 

which is obtained by the contraction of the symplectic 2-form u with the vector 
field X defining the fiow. 

For systems with finite, and furthermore even number of degrees of freedom, 
the symplectic 2-form is the inverse of the Hamiltonian structure functions which 
is the analog of the Hamiltonian operator. However, the generalization of the 
notion of an inverse to systems with infinitely many degrees of freedom is not 
immediate. That is, the symplectic 2-form is obtained by integrating the density 
|6| ^ 

a; = - du' A K^j du^ (75) 



18 



over the spatial variable, where K is the "inverse" of J. Given Hamiltonian 
operator J, its "inverse" may be defined by 

J'''Kkj=KjkJ''' ^5] (76) 

but eq.([76|) is an operator equation which acts on gradients of functions. On the 
other hand, quite generally there exist Casimir functions which are annihilated 
by the Hamiltonian operator. Such functions, c./. eq.(|37|) above, must be 



excluded in the definition of the "inverse" in eq.(76). 

Our approach to the construction of Hamiltonian operators is based on the 
construction of Dirac brackets for systems subject to second class primary con- 
straints. It is evident from eq. (p3|) that the essential element in the Dirac bracket 
is the inverse of the matrix of Poisson brackets of the constraints. On the other 
hand, according to eq. (^6|) we must invert again to obtain the symplectic 2-form, 
thus we have simply 

K,k{x)5{x ~ y) = {c^,{x)Ak{v)}, (77) 

that is, the inverse of the Hamiltonian operator can be obtained directly from 
the Poisson bracket of second class constraints. 



9.1 Symplectic 2-forms 



The Hamiltonian operator ( ^91 ) can be inverted in a straight-forward way subject 
to the provision above 



or we can immediately read it off from eqs.(|l7[) and (|77|). Eq. (psD is also a 
statement of the non-degeneracy of the Hamiltonian operator ( p9|) . Hence from 
eqs.([75|) and (|7|) we find the symplectic 2-form 

= — iPx dq - dp) A d{D-^p) (79) 

q 

corresponding to the local Hamiltonian operator (p9|). It can be verified that 
is a closed 2-form by direct calculation. We had found that in the u, q 
variables the inverse of Jq is a local operator and this is true for the family of 
local Hamiltonian operators ( p9|) in general. The symplectic 2-form assumes the 
simpler expression 



= \i^u^ (uxx dq - qx dux) A du (80) 



in these variables. By invoking the Poincare lemma, in a local neighborhood we 
can write 

= da^ = apx d{D^^p) (81) 
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and we note that the coefficient of the 1-form is also the Casimir ( |37| ) for the 
Hamiltonian operator J^. The closure of the symplectic 2-form ([7^) is equivalent 



to the satisfaction of the Jacobi identities by the Hamiltonian operator (29). 

For the nonlocal Hamiltonian operator ( |5l| ) the inverse is given by the local 
operator 

Qx^rr Qx^rr ^rr \ 

~jr \ " (82) 

Dx!^ J 

which also follows from eqs.(|47|) and (|77|). Then from eqs.([75|) and ( |8^ ) the 
symplectic 2-form appropriate to the nonlocal Hamiltonian operator is given by 

Px P^ J (83) 

= Krr dr A dp 

so that the verification that oj^ is closed is immediate. In a local neighborhood 
we can write 

uJk = da,^, = Kr dp (84) 

using the Poincare lemma. 



9.2 Symplectic form of equations of motion 

With the symplectic 2-forms ( [t^ ) and ( |8^ ) we need to check that eqs.(^ are 
satisfied. For this purpose we recall that given a 2-form uj — a{v,Vx) dv A dvx 
and the vector field X = m{v, Vx) d/dv, we have ixuj = (2amx + max) dv. If we 
consider luk with k — ^r^, from this expression we get 



Px 



Px 



Px 



Qx 

p7 



1 



] Po 

\ Px J 

_ 9x 

^Pxx 

: Px 



dp A dq^ 



dp - %dpx (85) 



Px 



where we have discarded a total derivative. Similarly, in order to check that 
ixijJii — dHfj^, we note that given the 2-form w = a{v,VxT--) dv A d(D~^v) 
and the vector field X — m{v,Vx, ■■■) d/dv, we have ix^^ = —amd{D,^^v) — 
mD~^{adv) . The application to the 2-form (|79|) yields 



ixi^i^i = -qxDx ^{(TqPxdq) + qxD^ ^(aqqxdp) 
= q<yq{Pxdq~ qxdp) = dHf_t 



(86) 



using eq.(pl). 
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9.3 Witten-Zuckerman 2-form 



Time plays a privileged role in Hamiltonian mechanics. While this presents no 
problem for systems with finitely many degrees of freedom, in field theory it has 
the disadvantage of non-covariance. In order to remedy this situation Crnkovic 
and Witten and Zuckerman js] have introduced the conserved current 2-form 
which provides an elegant covariant formulation of Hamiltonian structure. For 
Monge- Ampere covariance is particularly necessary because, as we noted in the 
introduction, the choice of time coordinate for RHMA2 is quite arbitrary. 

The simplest way to obtain the Witten-Zuckerman current 2-form u for the 
Lagrangian ( pH ) is to first construct the 1-form a which follows from the first 
variation of the Lagrangian and a is related to uj as in eg. ([7^) . Assuming the 
equations of motion m), the first variation reduces to a conservation law 



where 



SCx^a\, + a\,, (87) 

(88) 

Finally, the current 2-form ujx for RHMA2 is given by 

uj\ = Sa* = Xu^u^ {q^ Su A Su^ - u^^ Su A Sq) 

/ ,2 X (89) 

t^A = Sa"" Xu^u^ [-ijj^6u A + q^Su A 6q 



and using the equations of motion we can readily verify that it satisfies 

= 0, = 0, (90) 

where a ranges over two values t and x. The Witten-Zuckerman 2-form is closed 
and conserved. For the Lagrangian (^3|) a similar procedure yields 

W*^ = Krr Sr A Sp 

(91) 

— r Krr Sp A Sr 

which also satisfies eqs.(|90|). A useful relation in checking the conservation law 
for the 2-form ( ^l| ) is the dKdV, or the Riemann equation form of RHMA2 

rt=rr^ (92) 

which follows from the results of |l8| . 

We note that the time components of the Witten-Zuckerman 2-forms ( p9| ) 
and ( |9l| ) for RHMA2 that follow from the Lagrangians ( [III ) and (|4^) arc precisely 
the symplectic 2-forms (|o|) and (^3|) respectively. The use of the notation 6 for 
d follows M and is restricted to this section only. 
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10 Lax Pair for RHMA2 



We have seen that the recursion operators of section |^ satisfy the Lax equation 
( |62[ ) but these are not useful Lax pairs. We need to cast RHMA2 into the form 
of a zero-curvature condition 

f/t - 14 - [U, V]^0 (93) 

which is the basic element in the solution of of completely integrable systems 
using the inverse scattering transform. The zero-curvature condition for the 
SL{2, ii)-valued connection 1-form given by the pair 

U-( ^ 2 ] V-( ^^2 

(94) 

provides such a formulation of RHMA2. But in this case the potential has a 
quadratic dependence on the spectral parameter A which has so far not been 
considered for an application of inverse scattering techniques. This U, V pair 
is therefore not immediately amenable to treatment by the method of inverse 
scattering. 



11 Multi-Hamiltonian structure of RMA2 

The infinite classes of Hamiltonian operators we have obtained for RHMA2 
reduce to the compatible pair Jq and Ji of Hamiltonian operators (|30| ) and ( ^ 
when we consider the Hamiltonian structure of RMA2 with non-zero constant 
right hand side in eq. (||) . The corresponding pair of symplectic 2-forms are also 
unchanged and the only modification comes in the conserved quantities. The 
infinite sequence of conserved Hamiltonians for RMA2 are those which reduce 
to the RHMA2 Hamiltonians for /i = ^q^ and k — in the limit K ^ 0. 
Thus the basic Hamiltonian densities entering into the Lenard-Magri scheme 



JoSH^ = JiSH^ (95) 



are given by 



= l:q^P. + KD,-^p, ^--^(qJ + K) (96) 

Z Z Px 

and the infinite sequences of conserved Hamiltonians which are in involution 
with respect to Poisson brackets defined by Jo and Ji are modified for RMA2. 
For example, by the application of the recursion operator (|6^) to we get 



1 

2px 



Kp. 



2 



(97) 
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which, up to a divergence, is the same as the RHMA2 Hamiltonian density in 
the sequence ( |7T| ) for X ^ 0. Repeated apphcation of the recursion operator 
(|66|) yields 



Pt qx - Qt Px 



-D. 



{px 



(98) 



for the RMA2 hierarchy of equations. 

The elliptic case of eq.(^) is equivalent to the equation for minimal surfaces 
while the hyperbolic case corresponds to the Born-Infeld equation . Through 
the appropriate change of variables the rich multi-Hamiltonian structure of the 
Born-Infeld equation carries over into RMA2 which includes the Hamilto- 
nian operators Jq and Ji. 



12 Ur-RHMA2 

The local and nonlocal family of Hamiltonian operators for RHMA2 have scalar 
counterparts for the Ur-RHMA2 equation (||) which can be written as 

ut = — (99) 

Ux 

in the form of an evolution equation. With the definition 

2k 

this equation can be identified as 

Ut + UU., = 

which is the dispersionless KdV, or Riemann equation. Ur-RHMA2 admits 
infinitely many conserved quantities 

Hn = K (100) 

as well as infinitely many local Hamiltonian operators 

J„ = J^D,-^ (101) 

^a; '^xx '^x '^xx 

provided the various constants are related by 

n = 2(a + l), 2a(a + l)(2a-|- = -fc 

s o th at the equation of motion (|99| ) assumes the form of Hamilton's equations 
(103). Scalar Hamiltonian operators of this type were first considered by Vino- 
gradov . The family of nonlocal Hamiltonian operators 

= kpuiD~^ui (102) 
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which is due to Sokolov |£T[ is also appropriate to the Ur-RHMA2 equation. In 
this case the conserved Hamiltonian densities are also given by eq.(lOC) but now 
the constants are related by 

n = -2/3, 2/3(2/3 + l)fc/3 (/3 + l)fc 

and we find that Ur-RHMA2 is cast into Hamiltonian form with 

Ut = Jq 5H2a+2 = Jp SH^2/3 (103) 



which again results in an infinite set of Hamiltonian structures. In eqs.( |lOl| ) and 
(102) we have the gq-components of the RHMA2 matrix Hamiltonian operators 
(29[) and ( pi] ) respectively. 

The recursion operator obtained by the composition of these Hamiltonian 
operators 

= -^D,^ (104) 

Uxx Ux Uxx Ux 

can be factored as in eq.(|6^) with 

^ ^7^^^^ ^^^^^ 

Ux Uxx Ux 

resulting in a Sheftel'-type recursion operator for Ur-RHMA2. 



13 RHMA in arbitrary dimension 

In order to write RHMA„ as a system of nonlinear evolution equations it will 
be useful to introduce a compact notation. For this purpose we shall consider 
the determinants of the n — 1 x n — 1 matrices 



/ 



A'^- = (-1) 



fc+i 



det 



92 



Mil 

U21 



Ulk 



Uln-1 
U2n-1 



\ 



(106) 



i-ln-l / 



obtained by deleting the 0*'* row and k^^ column in the matrix of second deriva- 
tives. The latter is indicated by a hat over the omitted terms. In particular, for 
A; = we have the Monge- Ampere operator in n — 1 dimensions 

A = -A" ^ 

which is a statement of nondegeneracy of RHMA„. The system of evolution 
equations for RHMA„ is given by 



Ut 



qt = 



(107) 
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where qi — dq/dx^ and henceforth we shaU reserve the index i to range over n — 1 
independent variables while continuing the use of the summation convention over 
repeated indices. Thus the vector field 



X^4 + 1,.A^^ 
ou A oq 

defines the flow for RHMA in n dimensions. 



(108) 



Eqs.(107) are cast into the form of Hamilton's equations with the Hamilto- 
nian operator 

/ 



J = 



q 
li'A 



V 



qA' 
H'A' 



D 



(109) 



where ^ is an arbitrary differentiable function of q alone and prime denotes 
derivative with respect to the argument. The Hamiltonian function is given by 



(110) 



and once again, there exist infinitely many conserved quantities ( |110| ) and 
Hamiltonian operators (109) associated with RHMA in n dimensions. How- 
ever, this is not the full extent of the Hamiltonian structure of eqs.(107) as 
we have not considered nonlocal operators, or the dependence of fi on other 
variables. Concerning the latter point we note that 



(111) 



is also conserved for the system (107). Hence the number of independent vari- 
ables entering into the arbitrary function fi can be increased considerably with 
an attendant incre ase in the number of Hamiltonian operators which is already 
infinite in eq.(109). 



In o rder to present the symplectic structure of eqs.(107) we need the inverse 
of (109) and subject to the provisions of section ||, we find that it is again a 
local operator. Then from eq.(|75|) we get the symplectic 2-form 



w = — (AMm a dui + Adq A du) 
q 



(112) 



which can be directly verified to be a closed 2-form. In a local neighborhood we 
can write it as the exterior derivative of a 1-form a which is given by 

a^aAdu (113) 

where a is again related to /i through eq. (|38|) and a A is the Casimir for the 
Hamiltonian operator (109). 



25 



14 Geodesic flow for CHMA 



The Hamiltonian structure of the geodesic flow for CHMA is very similar to 
that of RHMA. Semmes |jl^ has introduced the notion of geodesies on A/", the 
space of smooth real-valued functions on / x M where J is a real interval. For 
F G Af{I X M) and {ddFy ^ the vector field 

defines geodesies on M. The holomorphie exterior derivative is denoted by d. 
Here as well as in the following, it will be understood that volume forms on M 
enclosed by square parantheses automatically carry the Hodge star operator so 
that the result is a 0-form. The discussion of the symplectic structure of CHMA 
by Semmes is based on the Kahler 2-form 

^ddF 

which is not the relevant obj ect t hat emerges from an examination of the Hamil- 



tonian structure of the fiow ( 114 ) . Our approach to the problem of the geodesic 



flow for CHMA will be in the framework of dynamical systems with infinitely 



many degrees of freedom and the resulting symplectic 2-form is given in eq.(118). 
The advantage of our approach lies in the direct proof it furnishes for the com- 
plete integrability of the geodesic fiow for CHMA. 

The geodesic flow for CHMA satisfies Hamilton's equations 

2)=^^l)=J^SH. (115) 

where X denotes the vector field (|114| ). The Hamiltonian density is given by 

n, = ^i [{ddFT] (iie) 

and the Hamiltonian operator is 
/ 



Re ■ 



[9g A {ddFr~^ nqdq A (99^)"-^ ] 

[f^'[{ddFr] ''^\ ^'[(aSF)"]' ^ ^ f,'[{ddFyf 

(117)- 



where n — fi{q) is an arbitrary differentiable function of its argument. 



The inverse of Hamiltonian operator (117), subject to the restrictions of 
section ^, is again a local operator which yields the symplectic 2-form 



aj^^ll{Re (dF A [dq A (ddF)"-^ A B] dF) + i [(99^)"] dq A dF } 



(118) 
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and for integrable complex structure w can be simplified by expressing the ex- 
terior derivative in terms of d, d. The statement of the symplectic structure of 



the geodesic flow for CHMA is given by eq.(|7j). The 2-form (118) is closed as 
one can show readily by direct calculation. However, it is more instructive to 
note that by invoking the Poincare lemma in a local neighborhood we can write 

= da, , Uc = -a{q) UddFY] dF (119) 
n 

where cr satisfies eq.(^) again. 

There are infinitely many symplectic 2-forms, compatible Hamiltonian op- 
erators and conserved Hamiltonians for the geodesic flow for CHMA. 



15 Conclusion 

We have considered the niulti-Hamiltonian structure of various real homoge- 
neous Monge- Ampere equations and found that quite generally they admit in- 
finitely many such structures. In particular for RHMA2 we have shown that 
there exist infinitely many Hamiltonian operators of both the local and non- 
local variety. The simplest Hamiltonian operator of the latter type leads to 
the Kac-Moody algebra of vector fields and functions on the unit circle. For 
Ur-RHMA2 we have the scalar version of the RHMA2 Hamiltonian operators. 
Finally, wc have shown that local Hamiltonian operators are generic to RHMA„. 
Thus the real homogeneous Monge- Ampere equation is a system with infinitely 
many Hamiltonian, or symplectic structures in the theory of integrable systems 
in arbitrary dimension. 
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